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Abstract 

In the context of the entangled B°B° state produced at the T(45) resonance, 
we consider a modification of the usual quantum-mechanical time evolution 
with a dissipative term, which contains only one parameter denoted by A and 
respects complete positivity. In this way a decoherence effect is introduced 
in the time evolution of the 2-particle B°B° state, which becomes stronger 
with increasing distance between the two particles. While our model of time 
evolution has decoherence for the 2-particle system, we assume that, after 
the decay of one of the two B mesons, the resulting 1-particle state obeys 
the purely quantum-mechanical time evolution. Prom the data on dilepton 
events we derive an upper bound on A. We also show how A is related to 
the so-called "decoherence parameter" £, which parameterizes decoherence in 
neutral flavoured meson-antimeson systems. 
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I. THE INTRODUCTION 



There is increasing interest in the last years in using particle physics phenomena for the 
study of possible deviations from quantum mechanics (QM). Efforts have been concentrated 
on two types of phenomena: oscillations, like K°-K° |l] and neutrino oscillations f2|, and 
quantum entanglement, where particularly suitable systems are the entangled K°K° and 
B°B° states || which are produced in e + e~ collisions at the resonances $ and T(4S"), 
respectively. These states become macroscopically extended objects before they decay. Thus 
in both types of phenomena macroscopic distances are involved. Furthermore, entangled 
systems are - due to EPR-Bell correlations || - important objects to clearly test QM against 
local realistic theories. Whereas entangled systems like K°K° have a rather longstanding 
and venerable position in the QM literature 0] , the physics of neutrino oscillations is a rather 
recent testing ground for QM; this development has been boosted, in particular, by the now 
well-established atmospheric neutrino anomaly, but also solar neutrinos and neutrinos in the 
early universe are discussed in this context. 

In this paper we concentrate on possible decoherence effects which might arise due to 
some fundamental modification of QM or due to the interaction of the system with its 
"environment" , whatever this may be. In the latter case, the idea of the influence of quantum 
gravity |||7J - quantum fluctuations in the space-time structure at the Planck mass scale - 
is especially attractive nowadays. Possible effects of the environment have been investigated 
intensively in the K°K° system in Refs. [|IJ^|-0]. But also other models of decoherence, like 
those found in Refs. [15|-p^|, may serve as working hypothesis. Our model, which we will 



propose in this paper, has some remote similarity with the models mentioned here, but ours 
will be tailored to the situation of two particles moving apart in their center of mass system. 

In the past, a measure of decoherence for entangled systems has been introduced on 
pure phenomenological grounds in order to determine quantitatively deviations from pure 
QM. This simple procedure of multiplying the quantum-mechanical interference term by 
1 — C 101 ' wnere C is called decoherence parameter, is a basis-dependent concept and works 
very well as a measure for interpolating continuously between pure QM (( = 0) and total 
decoherence (( = 1). The latter case corresponds to spontaneous factorization, also called 



Furry's hypothesis [pOf] . By investigating certain observables, the authors of Refs. [p21|-|25 



could show that the entangled K°K° and B°B° systems are far from total decoherence, at 
least when ( is introduced in relation to the basis of mass eigenstates, so that local realistic 
theories are highly unlikely. In other words, the presence of the interference term is well 



established in agreement with QM (see also Ref. [Q), which means that there is quantum 
interference of massive particles over macroscopic distances. 

In this paper we want to present a model of dissipation for entangled systems of two 
particles. In contrast to the prevailing concept in the literature, where dissipation is in- 
troduced at the 1-particle level and transferred to the 2-particle level through the tensor 
product structure of the Hilbert space of states (see, e.g., Ref. fl2fl), we assume the usual 



quantum-mechanical time evolution for the 1-particle states. Thereby we have in mind that 
entangled 2-particle systems become decoherent when they move apart over macroscopic 
distances, whereas for a 1-particle system QM is not modified. Our dissipative term in the 
2-particle time evolution obeys the condition of complete positivity P?fl . For reasons given 
below, we consider the entangled B°B° system with negative C parity. By using the exper- 
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imental value of the ratio R of the number of like-sign dilepton events over opposite-sign 
dilepton events, we can derive a bound on the strength of the dissipative term. Considering 
the time-integrated dilepton event rates, our model reproduces precisely the corresponding 
calculations with the phenomenological decoherence parameter ( associated with the Bh- 
Bl basis. As a result, we obtain a remarkably simple formula which relates the dissipative 
strength to the decoherence parameter (. In the context of the observable R, we also com- 
pare our model of 2-particle decoherence with the case where the analogous dissipative term 
is introduced already at the 1-particle level. 



II. THE MODEL 

Before considering the B°B° system, let us first discuss our model of decoherence in a 
2-dimensional Hilbert space of states Ti = C 2 . We allow for a non-hermitian Hamiltonian H, 
in order to include the possibility of incorporating particle decay in the Weisskopf-Wigner 



approximation f28|. We denote the normalized energy eigenstates by \ej) (j = 1,2) and 
have, therefore, 

H \ e j) = X j\ e j) with X j = m j ~ if] > (2- 1 ) 

where rrij and Tj are real and the latter quantities are positive in addition. Furthermore, 
we make the crucial assumption that 

(ei|e 2 ) = (2.2) 

despite the non-hermiticity of H. Including decoherence, the time evolution of the density 
matrix p has the form 

^ = -iHp + ipH^ - D[p\. (2.3) 

Our model of decoherence consists in assuming that 

D[p] = X (P lP P 2 + P 2 pP 1 ) , where P 3 = | ej )(e,| (2.4) 

and A is a positive constant. Such a term is also employed, for instance, in the context 
of neutrinos in the early universe (see, e.g., Ref. p9|). It can readily be checked that the 



decoherence term in Eq. ( |2.4| ) is of the Lindblad type |27J 



d\p\ = \ (E 4v + p E A ) A A - E a*p a ) > ( 2 - 5 ) 

V j 3 / 3 

if we make the identification A, = \f~\Pj. Thus, the term (|2~4l ) generates a completely 
positive map; moreover, since Pj = Pj and [Pj, H] = 0, the decoherence term would increase 
the "von Neumann entropy" and conserve energy in the case of a hermitian Hamiltonian 
(see Ref. |30| and references therein). However, what is more important in our discussion is 
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the fact that with the choice ( [2.4| ) the equations for the components of p decouple. Indeed, 
with 



P = Yl Pjk\ej)(e k \ , (2.6) 
j,k=i 

where pjk = p* k j, and with the time evolution Q2.3Q , we obtain 
Pn(*) = Pii(0)exp(-ri*), 

P22W =p 22 (0)exp(-r 2 t), (2.7) 
Pn{t) = Pi 2 (0) exp {~[i{m 1 - m 2 ) + (T 1 + T 2 )/2 + A] t} . 

Let us for a moment dwell on the motivation for our model of decoherence. To this end 



we start with the more general setting Aj = ■JXjPj (j = 1,2 and Xj > 0) with Pj = \pj){pj 
where the normalized vectors \pj) are linearly independent and in general different from the 
eigenvectors of the Hamiltonian. Note that we also allow for (px\p2) 7^ 0. In any case, one 
can use formula (|2.5|) to obtain a completely positive dissipative term in the time evolution 
( p.3|) , but in general one will not obtain the form of D[p] given by Eq. ( |2.4|) . For the time 
being we want to assume that H* = H 7^ holds and that H is non-degenerate in order to 
avoid trivial considerations. Now we have two possibilities: 1. The system {\px), |p 2 )} is the 
system of eigenvectors of H, i.e., it is equivalent to the orthonormal system {|ei), |e 2 )}; 2. 
{\Pi)i IP2)} is not equivalent to the system of eigenvectors of H. In the first case one can 
show that the form (|2.4|) of D[p] with A = (Ai + A 2 )/2 is obtained and that 



Case 1 & [H, Pj] = for j = 1, 2 : lim p(t) = Pip{0)Pi + P 2 p(0)P 2 (2.8) 



t—>oo 



holds. Furthermore, density matrices Pj, or linear combinations thereof, are constant solu- 
tions of the time evolution equation. In the second case, at least one of the vectors |p 2 ) 
is not an eigenvector of H and one can prove that 

Case 2 & 3 j = 1 or 2 with [H, PA ^ : lim p{t) = ±1 , (2.9) 

independent of p(0). We will see in the following - when we apply our model to the B°B° 
system - that the first case is closer to our physical intuition (see also last paragraph of 
this section). The two cases have been described in Refs. [^T],§2J in the context of 1-particle 
decoherence in neutrino oscillations. Note that Case 1 is used in Ref. (see also Ref. [ |3Cf 



whereas Case 2 is considered, e.g., in Ref. in the same context. If we allow for W ^ H, 
the picture, we have developed here, gets blurred because then there is a competition between 
particle decay, i.e., lim^oo Tr p(t) = 0, and the effect of decoherence. We nevertheless stick 
to the first case. Note that identifying the orthonormal system given by Pj (j = 1, 2) with 
the system of eigenvectors of H is not only motivated by the considerations above but also by 
simplicity; as we have seen in Eq. ( j2.7| ) we have decoupled time evolutions as a consequence. 
We want to stress, however, that in the case of CP violation, which is particularly important 
for the K°K° system, it might be useful to allow for small deviations from (ei|e 2 ) = 0, and 
thus for small deviations of the eigenvectors of the Hamiltonian from the orthonormal system 

{|Pl)> |P2>}- 
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Returning from general considerations, we now we apply our model of decoherence to 
the case of the 2-particle B°B° state, generated by the decay of the T(4S) resonance (for 
the formalism used in the B°B° system see, e.g., Ref. ||34||). We conceive t as the eigentime 
of B° and B° and make the identification 

\e l ) = \B H ®B L ) and |e 2 ) = \B L ® B H ) , (2.10) 

where the heavy and the light neutral B states are defined via 

\B H ) =p\B°) + q\B°) and \B L ) = p\B°) - q\B°) , (2.11) 

which have eigenvalues 

i % 
\h = m H - -T H and X L = m L - -T L , (2.12) 

respectively, of the effective 1-particle Hamiltonian Hi. For the 2-particle system, we trans- 
fer, as usual, the 1-particle Hamiltonian to the tensor product of the 1-particle Hilbert 
spaces by using H = Hi ® 1 + 1 ® ifi. We imagine that the first factor in the tensor product 
corresponds to particles moving to the left, whereas the second factor in the tensor prod- 
uct corresponds to right-moving particles. We assume CP conservation in B°-B° mixing, 
which is a good approximation p3)|35H and corresponds to \p/q\ = 1. In this case we have 
(B H \B L ) = and, therefore, (ei|e 2 ) = 0. In the following we will set p = q = 1/V2. At the 
T(4S') resonance, at t — 0, the entangled state 

M = ^ (K) - h» (2.13) 
is produced, which is equivalent to the density matrix 

P(0) = \ (|e 1 )(e 1 | + |e 2 )(e 2 | - | ei )(e 2 | - |e 2 )< ei |) . (2.14) 



With the time evolution ( |2.7| ), the initial condition (|2.14|) and taking into account that in 



the case of the vectors (|2.10|) we have Ai = A 2 = m H + m L — iY with Y = (r# + T L )/2, we 
obtain the time evolution 

P(t) = \^ 2Tt {|ei>(ei| + |e 2 )(e 2 | - e~ xt (| ei )(e 2 | + |e 2 )( ei |)} . (2.15) 



Note that the factor exp(— Xt) in the density matrix (|2.15| ) introduces decoherence as a 



consequence of the D-term in the time evolution ( p.3|) . In other words, for t > and A > 0, 
the density matrix ( |2.15| ) does not correspond to a pure state anymore. 

Having chosen the energy eigenstates (|2.10|) for the construction of the projectors Pj, our 
model complies with Case 1 ( |2.8| ). In this case, we would have no decoherence if p(0) = Pi 
or P 2 , though such initial conditions might be unrealistic. This agrees with our intention 
because in these cases we have no entanglement over macroscopic distances and no reason 
for modifying the quantum-mechanical time evolution. Note that using the projector states 
( p.lO| ) confines our Hilbert space of states to a 2-dimensional one. Using projector states 
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which are non-trivial orthogonal linear combinations of the states (|2.10 ), would lead to Case 
2 (|2.9p , still with a 2- dimensional Hilbert space. However, using projector states which are 
not linear combinations of the states ( |2.10| ), like, e.g., \B°®B°), \B°(B)B ), entails not only 
a time evolution into the full 4-dimensional Hilbert space of states, including \Bh ® Bu) 
and \Bl ® Bl), but also opens up the possibility for more involved schemes of decoherence 
than given by our simple model. 



III. THE MEASUREMENT 

In order to obtain information on the parameter A, which modifies the time evolution in 
the B°B° system, we adopt the following philosophy. We start at t = with the density 
matrix ( [2.14| ) for a B°B° state with negative C parity. This 2-particle density matrix 



follows the time evolution (2.15) and undergoes thereby some decoherence. We imagine a 



measurement of the B quantum number of the left-moving particle at time te and of the 
right-moving particle at time t r . For times min(te,t r ) <t< ma,x(te,t r ) we have a 1-particle 
state which we assume to evolve exactly according to QM, with the time evolution given by 
H x . 

In a mathematical language, we do the following. Assuming for definiteness tp < t r , at 
t = ti we calculate the trace 

Ti e {(\n)(n\®l)p(t e )}=p r (t e ;t = t e ), (3.1) 

where Tr^ means the trace evaluated only in the space of the left-moving particles; p(te) is 
given by Eq. Q2.15| ), evaluated at t — U; moreover, we have defined |1) = \B°) and |2) = \B°), 
and n = 1,2. Consequently, p r (tf, t = te) is a 1-particle density matrix for the right-moving 
ones. For t > te, it is denoted by p r (te;t) and follows the 1-particle time evolution. At 
t = t r , where we measure the B quantum number of the right-moving particles, we finally 
have (n' = 1, 2) 

N(n,tf,n',t r ) = Tr {\n')(n'\p r (tf, t r )} . (3.2) 
Using all the above formalism and allowing also for te > t r , we arrive at 

N(n,te;n',t r ) = ^e- r ^ 

x {\(n\B H )\ 2 \ (n'\B L )\ 2 e - Ar ^)/ 2 + | (n| B L ) | 2 | (n'\ B H ) \ 2 e Ar ^-^ 2 

_g-AniiBfe.tr) (( n \B H )(n\B L )*(n'\B L )(n'\B H )*e- iAm{ti - tr) 

+ (n\B L )(n\B H )*(n'\B H )(n'\B L )*e tAm ^)} . (3.3) 

In this equation we have used the notation Ar = Th — and Am = ran — mi- For the 
sake of clarity, we have retained the scalar products {ti\Bh,l) arid (n'\BH,L)- According to 
our assumption of CP conservation in B°-B° mixing, we will replace them by their values 
±l/y/2. It is easy to check that for A = one obtains the usual expressions found in the 
literature. Note that for te = t r and n = n' we have 

N(n, t e ; n, t f ) = V 2r ^ (l - e~^) , (3.4) 

which is different from zero, in contrast to the standard quantum-mechanical case. 
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IV. THE DILEPTONIC DECAYS 



In practice, measurement of the B quantum number of neutral mesons in the entangled 
B°B° state proceeds via flavour tagging when the mesons decay. Assuming the validity of 
the AB = AQ rule, in inclusive semileptonic decays £ + tags B° and i~ tags B° (£ = e 



or fx). In the following we will concentrate on dilepton events [5fJ]. Denoting the inclusive 
semileptonic decay rate by R», the numbers of dilepton events from the decay of ( |2.13 ) 
are then given by the integrals 



poo poo 

N ++ = Tj / dt t / dt r N(l,t t ',l,t r ), 
Jo Jo 

poo poo 

A__ = V] I dt t dt r N(2, ti ; 2, t r ) , 



JO JO 

poo poo 

N + _ = JV_+ = Tj dti \ dt r N(l, t e ; 2, t r ) . (4.1) 
Jo Jo 

Defining x = Am/T and y = Ar/2r and calculating these integrals leads to the result 

r| r i i 



N+ ~ = N ~+ = tr? \T^f + TT^ ^ ~ / • < 4 ' 3 ' 

where the function C(^) is given by the simple expression 

C(A) = with A = - . (4.4) 

sv j 2 + A r V ) 

It is interesting to note that x does not enter into £. 

Eqs. ( |4.2[) and ( |4.3|) reproduce the results of Refs. p2]j23[| , where the "decoherence pa- 



rameter" ( [|1^] is introduced phenomenologically in the observable 

R - N++ + N - (4 <i) 

by multiplying the interference terms in N ++ , etc., with (1 — £). In the model presented 
here, ( is expressed by the parameter A (see Eq. ( |4.4j )), the strength of the dissipative term 
in the modified time evolution ( |2.3j ). It has been discussed in the literature that the above 
phenomenological procedure of introducing a parameter ( depends on the basis chosen in 
B°-B° space p0|j22] - p4| , ^7i ■ In the present model it is the ( associated with the Bh-Bl basis. 



Let us perform a numerical estimate of £ and A along the lines presented in Ref. JzT 
To this end we use R (|4.5| ), which has been measured by the ARGUS and CLEO |39 
Collaborations. Combining both measurements, we obtain the value -R ex p — 0.189 ± 0.044 



[ 2jJ . As far as x is concerned, we use the value x exp = 0.740 ± 0.031 obtained by combining 



the data from all LEP experiments [0. With the approximation y = in Eqs. (|4.2j) and 



3|) |35],[40| and using the law of propagation of errors, from R exp and x exp we derive the 



following numerical estimates: 



7 



C = -0.06 ±0.10 and A = -0.11 ± 0.18. 



(4.6) 



The Belle Collaboration has published data on the correlated semileptonic decay rate as a 
function of the difference t(, — t r |41| . Of course, these data could also be used to put a limit 



on A, if we integrate N(n, if, n', t r ), Eq. (|3.3| ), over ti + t r . However, we do not have enough 
information to perform such a fit. 

Let us now compare our model, where decoherence is implemented at the 2-particle level, 
with the case where we have the analogous time evolution (K3) at the 1-particle level [T2| , |r3| ]. 
We use the same structure of the D-term as given by Eq. (|2.4j ), but now with |ei) = \Bu) 
and |e 2 ) = \B L ), instead of Eq. ( [2.1(J[ ). We denote the strength of the dissipative term by £, 



in order to distinguish it from A in the case of 2-particle decoherence. Evidently, we have the 
same time evolution fl2.7[) at the 1-particle level, with A replaced by £. Following the steps 
to derive Eq. (|3.3|), we obtain the same formula, except that exp(— A min(^, t r )) is replaced 
by exp(— £(te + t r )). Eventually, we arrive at N ++ and A + _ given by Eqs. (fL2|) and Q4.3|) , 
respectively, where ( is now given by[] 

c(s,x = A 7-1, 2 with E = ^- ( 4 - 7 ) 

Thus the two models of decoherence cannot be distinguished on the basis of the time- 
integrated dilepton event rates, but only on the basis the time-dependent event rates. A 
numerical estimate analogous to the one performed for A leads to the result H = — 0.04±0.07. 



V. THE SUMMARY 



In this paper we have considered a model of decoherence applicable in the center of mass 
system of two particles. Our model reflects the idea that, when the two particles move apart 
and eventually become macroscopically separated, some "forces" might be operative which 
de-entangle the quantum-mechanical state as a function of the distance. The dissipative 
term which we have added to the quantum-mechanical time evolution could be an effective 
term originating in some modification of QM; it could as well be based on some effective 
quantum-mechanical description of an interaction of the 2-particle system with an unknown 
environment. Our dissipative term respects complete positivity, which - we believe - is a 
useful physical guiding principle for modifications of the quantum-mechanical time evolution. 
In compliance with our idea, we assume that, after one of the particles has decayed, the 
other one follows the quantum-mechanical time evolution. We have applied our model in 
the case of the entangled B°B° state with negative C parity, where we have used the data 
on the B lifetime, the Bh-Bl mass difference measured by observing the time evolution of 
single neutral B mesons, and the ratio of like-sign over opposite-sign dilepton event rates 



x Note that the result for R with £ given by Eq. ( [4.71 ) agrees with the result for R in first order in 
H of Ref. [12 1, if the general Lindblad term in this paper is specialized to energy conservation. We 
thank the referee for pointing this out to us. 
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for the purpose of estimating the strength A of the dissipative term. In the case of time- 
integrated dilepton events, our simple model leads to a result which is also obtained by the 
phenomenological introduction of a "decoherence parameter" ( in the quantum-mechanical 
interference terms of the quantities N ++ , iV + _, etc. In the dissipative term D[p] we employ 
the states \Bh®Bl) and \Bl®Bh)] this, eventually, modifies the interference terms of N ++ , 



iV H , etc., with the £ associated with the Bh-Bl basis |p2|,p3| . Note that we have neglected 

CP violation in B°-B° mixing, which is a good approximation in this system. Transferring 
our model of decoherence to the K°K° system is not straightforward, because it requires to 
take CP violation and the non-orthogonality of the K s and K L states into account. Work 
on this is in progress. 
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